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Abstract 

' We study the spectroscopy and dominant decays of the bottomonium-like tetraquarks (bound diquarks-antidiquarks), focusing 

on the lowest lying P-wave [fo(/][bg] states Y\pq^ (with q = u,d), having J^'~^' = 1~~. To search for them, we analyse the BaBar 

' f _ 
data 1] obtained during an energy scan of the e'^e^ — )■ bb cross section in the range of -y/s = 10.54 to 11.20 GeV. We 



on the lowest lying P-wave [fo(/][bg] states yj;,,] (with q = u,d), having J^^ = 1 .To search for them, we analyse the BaBar 

Oh 

^ find that these data are consistent with the presence of an additional bb state Yjbq] with a mass of 10.90 GeV and a width 
of about 30 MeV apart from the T(5S') and T(65') resonances. A closeup of the energy region around the Yj^qj-mass may 
resolve this state in terms of the two mass eigenstates, ;] and V];, /!], with a mass difference, estimated as about 6 MeV. 

J — ■ We tentatively identify the state Yfj,,] (10900) from the _R6-scan with the state y6(10890) observed by Belle Q] in the process 

OO . e+e" ^ n(10890) ^ T(15', 2S) Tv^TT due to their proximity in masses and decay widths. 
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I. INTRODUCTION 



In the past several years, experiments at the two B-factories, BaBar and Belle, and at the Tevatron collider, CDF and 
DO, have discovered an impressive number of new hadronic states in the mass region of the charmonia These states 



generically labelled as X , Y and Z, however, defy a conventional cc charmonium interpretation 0, Moreover, they 
are quite numerous, with some 14 of them discovered by the last count, ranging in mass from the J^'~^' = 1++ X(3872), 
decaying into DD* , J/i/jtt^tt" , Jip'^, to the J^'~' = 1^^ F(4660), decaying into tA'tt+tt^ (for a recent experimental 
summary and references, see [6]). There is also evidence for an ss bound state, ys(2175) having the quantum numbers 
first observed by BaBar in the initial state radiation (ISR) process e+e~ — > 7isr /o(98O)0(1O2O), where 



/o(980) is the 0++ scalar state . This was later confirmed by BES S] and Belle 

These states are the subject of intense phenomenological studies. Three difi^erent frameworks have been suggested 
to accommodate them: (i) D — D* molecules 10-12|; (ii) ccg hybrids and (iii) Diquark-antidiquark or four quark 



states |14| - |16j . Of these hypotheses (i) and (iii) are more popular. For example, the motivation to explain the state 
X (3872), first observed by Belle Ll7]and later confirmed by CDF Q, DO Q and BaBar 1201 , as a hadronic molecule 
is that the mass of this state is very close to the D^D*^ threshold. Hence, in this picture, the binding energy is small 
implying that these are not compact hadrons, which have typical sizes of 0(1) Fermi. This makes it unlikely that such 
a loosely bound state could be produced promptly (i.e. not from B decays, as seen by Belle and BaBar) in high energy 
hadron collisions, unless one tailors the wave functions to avoid this conclusion. In particular, Bignamini et al. 21 1 
have estimated the prompt production cross section of X (3872) at the Tevatron, assuming it as a D^D*'^ hadron 
molecule. Their upper bound on the cross section pp — >■ X(3872) + ... is about two orders of magnitude smaller than 
the minimum production cross section from the CDF data 23|, disfavouring the molecular interpretation of X(3872). 



However, a dissenting estimate [22j] yields a much larger cross section, invoking the charm meson rescatterings. 

The case that the X, y, Z and are diquark-antidiquark hadrons, in which the diquark (antidiquark) pairs are in 
colour 3c (3c) configuration bound together by the QCD colour forces, has been forcefully made by Maiani, Polosa 
and their collaborators 14|-|l6j. The idea itself that diquarks in this colour configuration can play a fundamental 



role in hadron spectroscopy is rather old, going back well over thirty years to the suggestions by Jaffe 2j|. More 
recently, diquarks were revived by Jaffe and Wilczek [251] in the context of exotic hadron spectroscopy, in particu- 
lar, pentaquark baryons (antidiquark-antidiquark-quark), which now seem to have receded into oblivion. However, 
diquarks as constituents of hadronic matter may (eventually) find their rightful place in particle physics. Lately, in- 
terest in this proposal has re-emerged, with a well-founded theoretical interpretation of the low lying scalar mesons as 
dominantly diquark-antidiquark states and the ones lying higher in mass in the 1-2 GeV region as being dominantly 
qq mesons 26| . Evidence in favour of an attractive diquark (antidiquark) qq channel for the so-called good diquarks 



(colour antitriplet 3c, flavour antisymmetric 3/, spin-singlet positive parity) in the characterisation of Jaffe [4|] is now 



also emerging from more than one Lattice QCD studies 



27, 



for the light quark systems. On the other hand, no 



evidence is found on the lattice for an attractive diquark channel for the so-called bad diquarks (i.e., spin-1 states) 
involving light quarks 28]. However, as the effective QCD Lagrangian is spin-independent in the heavy quark limit, 
we anticipate that also the bad diquarks will be found to be in attractive channel for the [cq] and \bq] diquarks having 
a charm or a beauty quark. This implies a huge number of heavy tetraquark states, as we also show here for the 
hidden hh tetraquark spectroscopy. Earlier work along these lines has been reported in the literature using relativistic 



3CI. 



quark models [29| and QCD sum rules 

In this paper, we study the tetraquark picture in the bottom [bb) sector. In the first part (Section II), we classify 
these states according to their J^'-^ quantum numbers and calculate the mass spectrum of the diquarks-antidiquarks 
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[6(7] with q, q' = u, d, s and c in the ground and orbitally excited states by assuming both good and bad diquarks. 
The resulting mass spectrum for the 0++,l++,l^ , 1^^ and 2++ states having the valence diquark-antidiquark 
content [6(7] with q = u,d,s and c, and the mixed ones [5d][&s] (and charge conjugates) is shown in Fig. [TJ The 
main focus of this letter is on the J^'-' = 1^^ states, which are excited P-wave states. To be specific, there are four 
neutral states (n = 1, 4) with the quark content ([6u][6?2]) (which differ in their spin assignments) and another 
four ^[[^1 with the quark content ([M][6(J]). In the isospin symmetry limit, which is used in calculating the entries 
in Fig. [1] these mass states are degenerate for each n. Isospin-breaking introduces a mass splitting and the mass 
eigenstates called Yj["j and ^[["^^j (for lighter and heavier of the two) become linear combinations of Yj["j and Yj^^j"* . 
Thus, Yj["| = cos^? Yj["j + sin6' Yjl^^j and ^[1,"^] = ~ sin6' Yj[^j + cosO Y^j^y The mass differences are estimated to be 
small, with M(Yj^"'j) — Af(F[b ;]) = (7 ± 2) cos 20 MeV, where is a mixing angle. The electromagnetic couplings of 

(n) (n) 

the tetraquarks Yjj^ and Y^^ ^ are calculated assuming that the diquarks have point-like couplings with the photon, 
given by e(5[b<j] where e^/(47r) is the electromagnetic fine structure constant a and Qibq] — +1/3 for the [bu] and 
[6c] diquarks and Q[bq] = —2/3 for the [bd] and [6s] diquarks. Because of this charge assignment, electromagnetic 
couplings of the tetraquarks Yj["j and will depend on the mixing angle (Section III). 

To calculate the production cross sections e+e^ ^[t^i] ^ hadrons and e+e^ — > ^[i"'] ~^ hadrons, we need to 



calculate the partial widths rie^(Yj;, ;]) and ri"'' ( Y[f, jjj ) for decays into e+e^ pair and the hadronic decay widths 
r(Yj["j) and r(Fj[^"^^j). For the T{nS), the leptonic decay widths are determined by the wave functions at the origin 
^£'^5(0). The tetraquark states Yj^"j and ^j^"/]] are P-wave states, and we need the derivative of the corresponding wave 
functions at the origin, '^''^^^{0). To take into account the possibly larger hadronic size of the tetraquarks compared to 
that of the 66 mesons, we modify the Quarkonia potential, usually taken as a sum of linear (confining) and Coulombic 
(short-distance) parts. For example, the Buchmiiller-Tye QQ potential has the asymptotic forms V{r) ^ kqQ r 
(for r — > 00) and V{r) ^ 1/r ln(l/AQQj3 r^) (for r 0), where kqq is the string tension and Aqcd is the QCD scale 
parameter. The bound state tetraquark potential VgQ(r)^ will differ from the Quarkonia potential Vqqir) in the 
linear part, as the string tension in a diquark fegg is expected to be different than the corresponding string tension 
kqq in the QQ mesons, but as the diquarks- antidiquarks in the tetraquarks and the quarks-antiquarks in the mesons 
are in the same (3c3c) colour configuration, the Coulomb (short-distance) parts of the potentials will be similar. 
Defining k — kQQ/kqq, we expect k to have a value in the range k € [^, ^] 28 1. A value of k different from unity 
will modify the tetraquark wave functions ^^QQiO) from the corresponding ones of the bound 66 systems, effecting the 

(n) (n) 

leptonic decay widths of the tetraquarks. Hadronic decays of Yj), ;j and Y^f^ ^ are calculated by relating them to the 
corresponding decays of the T{5S), such as T(5S') — ?► B'^*^B'-*\ which we take from the PDG. We assume that the 
form factors in the two set of decays (Yjb,^] and T(5S')) are related by k, yielding the hadronic decay widths (Section 
IV). 

Having specified the mass spectrum and our dynamical assumptions for the tetraquark decays, we undertake a 
theoretical analysis of the existing data from BaBar [l| on Rb{s) = a{e^e^ bb)/a{e'^e^ /i+/i^), obtained during 
an energy scan of the e+e^ — > 66 cross section in the range of y/s = 10.54 to 11.20 GeV. The question that we ask and 
try to partially answer is: Are the kinematically allowed tetraquark states Y^'^^ and visible in the BaBar energy 
scan of Rb? To that end, we calculate the contributions of the lowest 1~~ tetraquark states Y[b,h] and Y^j ;] to the 
hadronic cross sections a^e'^e^ — > Y^b.i] — > hadrons) and (T(e^e^ — > Y^b.h] — > hadrons), and hence the corresponding 



^ We shall use the symbol Q and Q to denote a generic diquark and antidiquark, respectively. However, where the flavour content of the 
diquark is to be specified, we use the symbol [bq], and [bq] with q = u,d, s, c. 
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contributions Ai?f,(s).^ Our fits of tfie BaBar i?;,-data are consistent with tfie presence of a single state Yjbq] as a 
Breit-Wigner resonance with the mass around 10.90 GeV and a width of about 30 MeV, in addition to the T{5S) and 
T{6S). The quahty of the fit with three Breit-Wigners is found to be better than the one obtained with just 2 (i.e., 
with T{5S) and T{6S)), as reported by BaBar [l| (Section V). A closeup of the energy region around 10.90 GeV is 
necessary to confirm and resolve the structure reported by us, as the isospin-induced mass difference between the two 
eigenstates yj^.^j and Y[;, ;] comes out as about 6 MeV, which is comparable to the BaBar centre-of-mass energy step 
of 5 MeV. We hope that this can be investigated in the near future by Belle. 

We tentatively identify the state Y|fcq] (10900) with the state yf,(10890) measured in the process e+e^ 1^(10890) ^• 
T(1S', 25*) 7r+7r~ Q]. An analysis 32] of the Belle data on the decay widths r{Yb — > T{1S, 2S) tt+tt^), dipion invariant 
mass spectra and the helicity angular distributions is in agreement with the tetraquark interpretation presented here. 



II. SPECTRUM OF BOTTOM DIQUARK-ANTIDIQUARK STATES 

The mass spectrum of tetraquarks [bq] [bq'j with q = u, d, s and c can be described in terms of the constituent 
diquark masses, rriQ, spin-spin interactions inside the single diquark, spin-spin interaction between quark and antiquark 
belonging to two diquarks, spin-orbit, and purely orbital term 0, i.e. 

H^2mQ + H^!^/^ + Hfs^^ + Hsl + Hll, (1) 



where: 



41^^ = 2(/Cfc,)3[(S,.S,) + (S-,.S,)], 

-ffgf^' = 2(/Cfc5)(S6 • Sg + • Sg) + 2/Ch6(Sh • S^) + 2/C55(Sq • Sq), 

Hsl = 2Aq(Sq-L + S2-L), 

Hll = Sg^rf|S±l). (2) 

Here mg is the mass of the diquark [fog], (/Cfcq)^ is the spin-spin interaction between the quarks inside the diquarks, 
ICbq are the couplings ranging outside the diquark shells, Aq is the spin-orbit coupling of diquark and Bq corresponds 
to the contribution of the total angular momentum of the diquark- antidiquark system to its mass. The overall factor 
of 2 is used customarily in the literature. For the calculation of the masses we assume isospin symmetry, i.e. the 
isodoublet consisting of the states 

>li:J = H[H and yg = [6d][H (3) 

are degenerate in mass for each n. Later, we will calculate the isospin symmetry breaking effects in the masses. 

The parameters involved in the above Hamiltonian ([2]) can be obtained from the known meson and baryon masses 
by resorting to the constituent quark model jssi 

i7 = ^m, + ^2/C,, (S, -S,), (4) 

i i<j 

where the sum runs over the hadron constituents. The coefficient /Cy depends on the flavour of the constituents i, 
j and on the particular colour state of the pair. Using the entries in the PDG for hadron masses along with the 



^ We shall often refer to the ground states ^j^^^j and Yy^q without the superscript for ease of writing. 
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TABLE I: Constituent quark masses derived from the L — mesons and baryons. 



Constituent mass (MeV) 


q 


s 


c 


b 


Mesons 


305 


490 


1670 


5008 


Baryons 


362 


546 


1721 


5050 



TABLE IL Spin-Spin couplings for quark-antiquark pairs in the colour singlet state from the known mesons. 



Spin-spin couplings 


qq 


sq 


ss 


cq 


cs 


cc 


bq 


bs 


be 


bb 


(/C,,)„(MeV) 


318 


200 


129 


71 


72 


59 


23 


23 


20 


36 



assumption that the spin-spin interactions are independent of whether the quarks belong to a meson or a diquark, the 
results for diquark masses corresponding to X (3872) and Y (2175) were calculated in the literature Here, we 

extend this procedure to the tetraquarks [bq] \bq] . The constituent quark masses and the couplings /C^ for the colour 
singlet and antitriplet states are given in Table I, II and III. 

To calculate the spin-spin interaction of the QQ states explicitly, we use the non-relativistic notation |S'g, 5*^; J), 
where Sq and Sq are the spin of diquark and antidiquark, respectively, and J is the total angular momentum. These 
states are then defined in terms of the direct product of the 2x2 matrices in spinor space, T", which can be written 
in terms of the Pauli matrices as: 



0-2 



pO " ^ . -pi 

V2' 



1 



V2 



(5) 



which then lead to the following definitions: 



|0q,02; 


Oj) = 




Oj) = 




l./> - 








1./) - 



1 



{(72) (^2) , 



1 



((72(7') (g) ((72(7' 



^ ((72) <S) ((72(7 



((72(7*) <S) ((72) 



(72(7 



(6) 



The properties of these matrices are given in the appendix of ref. [l^ . The next step is the diagonalization of 
the Hamiltonian ([T]) using the basis of states with definite diquark and antidiquark spin and total angular momen- 
tum. There are two different possibilities [l^: Lowest lying [6(j'][6q] states {Lqq ~ 0) and higher mass states 
{Lqq = 1), which we discuss below. 



TABLE III: Spin- Spin couplings for quark-quark pairs in colour 3 state from the known baryons. 



Spin-Spin couplings 


qq 


sq 


cq 


cs 


ss 


bq 


bs 


be 


(/C,,)3(MeV) 


98 


65 


22 


24 


72 


6 


25 


10 
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A. Lowest lying states {Lqq = 0) 



The states can be classified in terms of the diquark and antidiquark spin, Sq and Sq, total angular momentum 
J, parity, P and charge conjugation, C. Considering both 1700^ and bad diquarks and having Lqq = we have six 
possible states which are listed below. 



i. Two states with J^"^ = 



|0++) = |OQ,Og; Oj); 

|0++') = 0,;). (7) 



ii. Three states with J = 1: 



= TfdOS'lg; l./> + |ls,Og; 1,/)); 
= ^{\OqAq, lj)-|lQ,Og; 1,/)); 

= |1q>1q; I./)- (8) 

All these states have positive parity as both the good and bad diquarks have positive parity and Lqq = 0. The 
difference is in the charge conjugation quantum number, the state is even under charge conjugation, whereas 

and are odd. 

iii. One state with J^'^ = 2++: 

|2++) = |iQ,ie; 2j). (9) 

Keeping in view that for Lqq — there is no spin-orbit and purely orbital term, the Hamiltonian ([l} takes the 
form 

+2(/C6,-)(Sb • S,- + S5 • S,) + 2/C,5(S6 • S5). (10) 

The diagonalisation of the Hamiltonian pUj) with the states defined above gives the eigenvalues which are needed to 
estimate the masses of these states. It is straightforward to see that for the 1++ and 2++ states the Hamiltonian is 
diagonal with the eigenvalues [l^l 

M (1++) = 2m[b,] - (/Cb,)3 + i/C,,- - ICbg + i/C,fc, (11) 
M{2++) = 2TO[,,] + (/Ct,)3 + i/C5,- + /Cb,-+i/C,6. (12) 

All other quantities are now specified except the mass of the constituent diquark. We take the Belle data Q as input 
and identify the y;,(10890) with the lightest of the 1^^ states, ^[6^], yielding a diquark mass TO[(,g] = 5.251 GeV. 
This procedure is analogous to what was done in pL4i], in which the mass of the diquark [cq] was fixed by using the 
mass of X{3872) as input, yielding m[j,q] = 1.933 GeV. Instead, if we use this determination of m[j,q] and use the 
formula TO[f,g] = + (to6 — m-c), which has the virtue that the mass difference nic — rrib is well determined, we get 
TO[b^] = 5.267 GeV, yielding a difference of 16 MeV. This can be taken as an estimate of the theoretical error on m^j^j , 
which then yields an uncertainty of about 30 MeV in the estimates of the tetraquark masses of interest for us. 
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The couplings corresponding to the spin-spin interactions have been calculated for the colour singlet and colour 
antitriplet only. In Eq. ([2]), however, the quantities /Cqq, ICbq and JCfj/j involve both colour singlet and colour octet 
couplings between the quarks and antiquraks in a QQ system. So for /C^j [3] 



- 1 2 



(13) 



where (/Cj{,)q is reported in Table II. (/Cjf,)g can be derived from the one gluon exchange model by using the relation 



with (X) = 0, 4/3, 4/3, 3 for X = 0, 3, 3, 8 respectively. Finally, Eq. ^ gives 

^hb (NN) = \ (^6b)o ■ 



(14) 



(15) 



Now, we have all the input parameters to calculate the mass spectrum numerically. Putting everything together 
the masses for the hidden bb tetraquark states 1++ and 2++ states are: 



M (1++) = 10.504 GeV, for q ^ u, d, 

= 10.849 GeV, for g = s, 

= 13.217 GeV, for q ^ c, 

M (2++) = 10.520 GeV, for q = u, d, 

^ 10.901 GeV, for q = s, 

= 13.239 GeV, for q ^ c. 



(16) 
(17) 
(18) 
(19) 
(20) 
(21) 



For the corresponding 0"*"^ and 1+ tetraquark states, the Hamiltonian is not diagonal and we have the following 
2x2 matrices: 



-3(/C, 



bq)3 



V3 



(22) 



Mil 



+-\ - 



^{^bq)^ + ^bq 



2 



ICqq — JCbb 



ICqq — K-bb 



{JC 



bq)3 



ICbq 



(23) 



To estimate the masses of these two states, one has to diagonalise the above matrices. After doing this, the mass 
spectrum of these hh states is shown in Fig. [T] 

B. Higher mass [feg][&g] states {Lqq — 1) 



We now discuss orbital excitations with Lqq — 1 having both good and bad diquarks. In this paper, we are 
particularly interested in the 1 multiplet. Using the basis vectors defined in reference [3l the mass shift due to 
the spin-spin interaction terms Hss becomes: 



/ —3 (/Cfjg), 







AMss = 





(fCbq)^ — JCbq + {ICqq + JCiji) /2 

~ (ICbq)^- ICbq - {ICqq + ICbb) /2 J 



(24) 
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TABLE IV: Eigenvalues of the spin-orbit and angular momentum operator in Eq. ([T} for the states having J — Lqq + Sqq — 1. 



I'S's 1 'S'q, Sqq, Lqq) 


a {Sq , Sq, Sqq, Lqq) 




|0, 0, 0, 1) 





1 


|1, 0, 1, 1) 


-2 


1 


|1, 1, 2, 1) 


-6 


1 


|1, 1, 1, 1) 


-2 


1 


|1, 1, 0, 1) 





1 



The eigenvalues of the spin-orbit and angular momentum operators given in Eq. ([!]) were calculated by Polosa et 
al. [l^l, and we have summarised these values in Table ITVl '^ 

Hence the eight tetraquark states [6q][6q] {q — u,d) having the quantum numbers 1 are: 



M^y' , {Sq = 0, 5g = 0, Sqq = 0, Lqq = 1) = 2m[bg] + Al + Bq, 



r(4) 



1, Sq 
1, Sq 



0, 5, 
1: S, 



QQ 



QQ 



1, L 
0, L 



QQ 



QQ 



2™ 



2m 



(^S - 1' 5q = 1' 2, Lqq = 1) = 2™[5,j 



A + As - 2ylg 4 

2A + A3 + BQ, 

2A + A3 - 6Aq 



Be 



Be 



(25) 



where Xi{i — 1,2,3) are the diagonal elements of the matrix AM55 given in Eq. ((24)) . Note that there are 16 
electrically neutral self-conjugate 1 tetraquark states Y^^i^^ with the quark contents with q = u,d,s or c, 

of which the two corresponding to and i.e., ^[Ii"! and Y^j^^^ are degenerate in mass due to the isospin 

symmetry. There are yet more electrically neutral J^^ = 1 states with the mixed light quark content [fed] [5s] and 
their charge conjugates [6s] [fed]. However, these mixed states don't couple directly to the photons, or the gluon, 
and are not of immediate interest to us in this paper. 

The numerical values of the coefficients corresponding to Aq and Bq are given in Table |IV] and are labelled by a 
and b, respectively. The quantity A is the mass difference of the good and the bad diquarks, i.e. 



A = mg {Sq = 1) - mg {Sq = 0) . 



(26) 



In order to calculate the numerical values of these states, we have to estimate A which is the only unknown remaining 
in this calculation. Following Jaffe and Wilczek the value of A for diquark [bq] is A = 202 MeV for q — u, d, s 
and c quarks. We recall that we have used the known mesons and baryons to calculate the couplings of the spin-spin 
interaction and we can extend the same procedure to the S* = 1, L = (0, 1) meson states B* , Bi (5721), B2 (5747) to 
calculate the values of Aq and Bq which are: 



Aq ^ 5 MeV, for q = u, d, 

Aq ^ ?, MeV, for q = s,c, 

Bq = 408 MeV, for q = u, d, 

Bq = 423 MeV, for q^ s,c. 



(27) 



^ The entry for a in the last row of Table ITVl differs from the corresponding one in the first reference in which is given as —2, but 
this point has now been settled amicably in favour of the value given here. 
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TABLE V: Masses of the 1 tetraquark states M^^j^ in GeV as computed from Eqs. (I25|, and (HTjl. The value 
(for q = u,d) is fixed to be 10.890 GeV, identifying this with the mass of the It from Belle Q| 



Ibq] 





q — u, d 


q = s 


g = c 


q — d, q — s 




10.890 


11.218 


13.618 


11.054 














11.130 


11.479 


13.841 


11.281 














11.257 


11.646 


14.025 


11.476 














11.227 


11.629 


14.009 


11.453 



Numerical values of the masses for the states given in Eq. psp are quot ed in Table V. Some of the entries, in particular 



(q = u, d, s), are comparable with the existing ones in refs. 
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[bl] 

Finally, the mass spectrum for the tetraquark states [fo^lifo?] for q = u,d,s,c with J^*-^ — 0++, 1^^, 1^ , 1^^ and 
2"'""'" states is plotted in Fig. [T]in the isospin-symmetry limit. The bb tetraquark states with mixed light quark content 
[6(i][6s] are also shown in this figure. Of these, the state Yjl^^j (10.890) shown in the upper left frame in Fig. [T]is 
of central interest to us in this paper. 

III. ISOSPIN BREAKING AND LEPTONIC DECAY WIDTHS OF THE J^^ = 1 " TETRAQUARKS 

We discuss in this section the isospin breaking effects, which were neglected in the previous section, and calculate 
the decay widths Tee{Y^b,i]) and Tee{Y^b.h]) for ^[6,;] and Y^i, ^. The mass eigenstates are given by a linear superposition 
of the states defined in ([3]). Introducing a mixing angle 9, we have, for the lighter and heavier states: 



Y, 



[b.i] 



cos Y[b„] + sin Y[M] , 



Y[b,h] = - sin 6* Y^hu] + cos 9 Fm] . 

The isospin breaking part of the mass matrix is 

/ 2m„ + (5 S \ 
\ S 2md + S I ' 



(28) 
(29) 



(30) 



where S is the contribution from quark annihilation diagrams, where the light quark pair annihilates to intermediate 
gluons. Taking this into account, the isospin mass breaking is given by** 



M{Y[bM]) - M{Y[bj]) - (7± 3) cos(2^) MeV. 



(31) 



The partial electronic widths TeeiY^bfi) and TedY^b^h]) are given by the well-known Van Royen-Weisskopf formula 
for the P-states, which we write generically as: 



167rQ2c,2|vI;'g_(0)|2 



(32) 



The expression I I31II differs from the one derived in [l^ . but there is consensus now on the expression given here. 
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FIG. 1: Tetraquark mass spectrum with the valence quark content with q — u,d, assuming isospin symmetry (upper left 

frame), with q = s (upper right frame), with q — c (lower left frame), and for the mixed light quark content [6d][&s] (lower right 
frame). Some important decay thresholds are indicated by dashed lines. The value 10890 is an input for the lowest J^^ = 1 
tetraquaxk state Yj^^j . All masses are given in MeV. 



where Q = Qibd] = —2/3 is the diquark charge in Ym — [bd]\bd] and Q = Q[bu] = +1/3 is the charge of the diquarks 
in Ybu = [bu\[bu\, and ^gg('^) = is the first derivative in r of the wave function of the tetraquark, which 

needs to be taken at the origin, i.e. ^Qg(O) = ■\/3/ (47r)i?'(0). We have approximated u by the diquark mass. 

We determine the wave functions for the P-state tetraquarks [bd] [bd] and [bu] [bu] from the corresponding wave 
functions for the P-state bb system by scaUng the string tension in the linear part of the potential, as discussed in the 
introduction. As most potential models agree in the., (c„„E„i.g, p„(s a.d the l..e„ „1 U,e p„«„«al 

essentially determines the heavy Quarkonia wave functions, the uncertainty in ^'^^(O) from the underlying model is 
not a concern. We have used the QQ-onia package of 3J], yielding |i?'(0)p = 2.062 GeV^ for the 66 radial wave 
function, which we have used as normalisation. The corresponding value for the tetraquark states [6(7] [6g] is then 
calculated as ^qq(O) — K^'f,b(0), and used in our derivations of the decay widths. We expect that for all the P-states 
Y^bu] ^^'^ ^[bd] ' ^^"^ electronic widths will be constant, to a good approximation. 
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The ratio Tlee{Yb) of Tee (^[f,,;]) and Tee{Yibji]) is given by 



7^ee(n; 



TeciYlb.h]) QUO) 



1 - 2tan6l 
2 + tan 6* 



(33) 



where Qi{0) — Q[bu] cos 6* + Q[bd] sinO and Qh{0) — ~Q[bu] sin 6* + Q[bd] cosO are the mixing-angle weighted charges. 
Since the total cross sections for e^e" — >■ (^jb./j , ^jb./i]) — > hadrons are directly proportional to Tee{Y[b.i]) and ree(Yjt, ^jj), 
the ratio TZee{Yb) is accessible from the experiment. The absolute values of the decay widths Tee{Y[b^i]) and TgeiY^b.h]) 
are given by TeeiY^bj]) = 0.4 K^Qi{6)^ keV, where Qi{0) are the mixing angle weighted charges of the two mass 
eigenstates, Y[b_i] and Y^b^h]- which can also be seen in ((33)) . 



IV. DIQUARK-ANTIDIQUARK DECAY MODES 

In this section we discuss the dominant hadronic decays of the Lqq — 1 states. In doing this, we restrict ourselves 
I n the two-body decays, Yjbgj -> Bq'Bq, and when allowed kinematically, also the decay Y^bq] ~> Af,Af,. Their 
thresholds are pictured in figure [1] These decays are Zweig allowed and involve essentially quark rearrangements and 
the possible pop-up of a light qq pair to make the AbAb state. The decays Y^bq] ^ T(15, 25) tt+tt" are also Zweig 
allowed. However, they are sub-dominant and can be neglected in estimating the total decay widths. 

The vertices and the corresponding decay widths of the dominant decays are given below: 




M P <^ 



F{g>^P{q + ir 
-gf'^k + q)P 
+gp''{q + kY) 



a/3 



r 



r 



'211.13 



F^\k 



F2|fc|3 
4M27r ■ 



(34) 



F^|A;|3(48|fcp-104M^|fcp-|-27M"^) 
27r(M3-4|fc|2M)2 ■ 



47r M^n 



The centre-of-mass momentum |fc| is given by 



\k\ 



^NP - (Ml + M2) V^'^^ - i^h -M2f 



2M 



(35) 
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where M is the mass of the decaying particle and Mi, M2 are the masses of the decay products. The matrix elements 
are obtained by multiplying the vertices in (p4)) by the polarisation vectors. Thus, for the decay Y[{,^g] — > BqBq, the 
Lorentz-invariant matrix element is given by = ej}'''''^ F{k^ — l^), and likewise for the other decays shown above. 
The decay constants F and F' are non-perturbative quantities. We estimate them using the known two-body decays 
of T{5S), which are described by the same vertices as given above. The different hadronic sizes of the bb Onia states 
and the tetraquarks Yj^^] are taken into account by the quantity k, discussed earlier. We use the partial decay widths 
for the decays T(5S') BB,BB*,B*B* from the PDG values of the fuh width, given as rtot[T(55)] = 110 ± 13 
MeV [3] and the respective branching ratios. They are called Fpdg and given in Table I VIII yielding the coupling 
constants, called i^pDGi and |fc|. For the decays Yj^^^j — > Af,Af, and Yjl^^^j — >■ 5S, we take F ^ F' = l.l^Qgg, and include 
a factor of 1/3 for the baryonic final state to take into account the creation of the qq pair from the vacuum. We 
remark that the estimates of Fpdg will be modified, if as anticipated by the BaBar i?f,-analysis P], the decay width 
rtot[f (55*)] has a significantly lower value. 

The input values for the masses used in our calculation are listed in Table IVTl With this input, our estimates of 
the decay widths for Yj^^^^'j are given in Table IVllIl We also give the total decay widths (up to the factor k^). As seen 
in this table, the lowest lying 1^^ states Yj[^| are expected to have decay widths of 0(50) MeV, for — 0.5. Thus, 
the decay widths of Yj^^j are consistent with the corresponding measurements by Belle, if we identify Yj^^^j with their 
Yf,. The higher 1 states have much larger decay widths and will be correspondingly more difhcult to find. 

TABLE VI: Input masses taken from Q| in units of GeV. 



hadron 


mass 


hadron 


mass 


hadron 


mass 


B 


5.279 


TT 


0.139 


T(IS') 


9.46 


B* 


5.325 


At, 


5.62 


T(4S) 


10.5794 


Bs 


5.366 


Hi, 


5.792 


T(10860) 


10.865 


Bt 


5.412 


K 


0.4937 


T(11020) 


11.019 



TABLE VII: 2-body decays T(5S') B^*>B^*>, which we use as a reference, with the mass and the decay widths taken from 
Q|. The extracted values of the coupling constants Fpdg and the centre of mass momentum are also shown. 



process 


rpDG[MeV] 


FpoG 


|fc|[GeV] 


T(10860) B B 


< 13.2 


< 2.15 


1.3 


T(10860) B B* 


i5.4«;6 


q 7+0.7 
■5- ' -0.9 


1.2 


T(10860) B* B* 


48t}l 


1 +0.13 
J^-0.12 


1.0 



V. ANALYSIS OF THE BABAR Rt ENERGY SCAN AND POSSIBLE SIGNAL OF A bb TETRAQUARK 
STATE AT 10.90 GEV 

BaBar has recently reported the e+e~ — >■ bb cross section measured in a dedicated energy scan in the range 10.54 
GeV and 11.20 GeV taken in steps of 5 MeV Their measurements are shown in Fig. [2] (left frame) together with the 
result of the BaBar fit, the details of which are described in their paper and which were also made available to us [s^. 
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TABLE VIII: Reduced partial decay widths for the tetraquarks ^[j^p the extracted value of the coupling constant F and the 



centre of mass momentum (top left). The reduced total decay widths for yi*', are also tabulated (top right) and for the 



tetraquarks vi'^', (the lower two tables). The errors in the entries correspond to the errors in the decay widths in Table IViT 



[6<7l 



Decay Mode 


r/«;2[MeV] 


F 


|fcj[GeV] 


^fM B 


< 15 


2.15 


1.3 






3.7 


1.2 


y,w ^ B* B* 




1 


1.1 


Yl'l ^B B 


< 33 


2.15 


1.8 


^fM ^B B* 


40+18 


3.7 


1.7 


^fU ^ B* B' 


I62I42 


1 


1.6 


^'[M ^B B 


< 43 


2.15 


2 


y(^j ^B B' 




3.7 


1.9 


^ B* 


Oqi +60 


1 


1.8 




lOl^ 


i-itSJl/s 


0.3 


^BB 


< 41 


2.15 


1.9 


y^^'^B B' 


541^3 


3.7 


1.8 


yi%^B' B- 


2131^^ 


1 


1.8 



1 Tetraquark 


rtot/«:2[MeV] 




88± 16 


v{2) 


238 ± 48 




342 ± 65 




308 ± 60 



Decay Mode 


r/«:2[MeV] 


F 


IfcKGeV] 


y^rj ^ Bs 


< 26 


2.15 


1.6 


^ll'} ^ Bs Bt 


00+14 

■JO_j4 


3.7 


1.6 


^ b: b: 


1181^° 


1 


1.5 




< 47 


2.15 


2 


^ Bs Bt 


64lf, 


3.7 


2 


-> -b: b: 


258t«^ 


1 


1.9 


y\bs] Bs Bs 


< 63 


2.15 


2.3 


^[1!] ^ b: 


86lf, 


3.7 


2.2 


^ b: 


3671^° 


1 


2.1 




1911:; 


1-11^35/3 


0.6 


yjltj -> Bs 


< 61 


2.15 


2.2 


^IM ^ -B. B: 


8411 


3.7 


2.2 


^[S ^ b: 


355lgQ 


1 


2.1 


'■[bs] ^ - - 


J^"-io 


^•^-0.35/ ■J 


0.5 



1 Tetraquark 


Ttot/n[MeV] 


^\bs] 


176 ± 33 


^\bs] 


368 ± 70 


^[bs] 


534 ± 100 


y-{4) 


516 ±96 



Their fit model of the i?fc-data contains the following ingredients: a flat component representing the 66-continuum 
states not interfering with resonant decays, called Anr, added incoherently to a second flat component, called A^, 
interfering with two relativistic Breit-Wigner resonances, having the amplitudes Aioseoi ^11020 and strong phases, 
(/>10860 and (/>iio20, respectively. Thus, 

a{e+e- ^ hi) = |A„,|2 + + Aio86oe*'^^"«™BVF(Mio860,rio86o) + ^ii02oe^'^"°^°BW^(Mno2o,riio2o)|' , (36) 
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with BW{M,T) ~ — M'^) + iAIT]. The results summarised in their Table II for the masses and widths of the 
T{5S) and T(6S') differ substantially from the corresponding PDG values in particular, for the widths, which are 
found to be 43±4 MeV for the T(10860), as against the PDG value of 110±13 MeV, and 37±2 MeV for the T(11020), 
as compared to 79 ± 16 MeV in PDG. As the systematic errors from the various thresholds are not taken into account, 
this mismatch needs further study. The fit shown in Fig. [2] (left frame) is not particularly impressive having a x^/d.o.f. 
of approximately 2. In particular, the data points around 10.89 GeV and 11.2 GeV lie systematically above the fit. 
In our analysis of the BaBar data, we were able to reproduce these features, but also found that the fit-quality can 
be improved somewhat at the expense of strong phases 0io86O and 0iio2O: which come out different than the ones 
reported by BaBar P|. We do not show this fit here as the resulting i?f,-line-shape is close to the one shown in the 
BaBar publication and reproduced here. 

We have repeated the fits of the BaBar i?b-data, modifying the fit model in Eq. (p6|) by taking into account two 
additional resonances, corresponding to the masses and widths of ^[f,.;] and Y^i, ^. Thus, formula (1361) is extended by 
two more terms 

^y,,,„e^^-[MiW(A^y,, „,ryj, ,j) and A^.,,,,, e^^-i'.'-i SVF(Afy,,,„ , Fyj. ^j ), (37) 
which interfere with the resonant amplitude Ar and the two resonant amplitudes for T(55') and T(6S') shown in 

n 

Eq. (|36p . We use the same non-resonant amplitude Anr and Ar as in the BaBar analysis [1]. The resulting fit 
is shown in Fig. [2] (right frame). Values of the best-fit parameters are shown in Table IIXI from where one see 
that the masses of the T{5S) and T(6S') and their respective full widths from our fit are almost identical to the 
values obtained by BaBar [l|. However, quite strikingly, a third resonances is seen in the i?b-line-shape at a mass 
of 10.90 GeV, tantalisingly close to the lb(10890)-mass in the Belle measurement of the cross section for e~^e~' — 
yb(10890) T(15, 25) tt+tt-, and a width of about 28 MeV. In the region around 11.15 GeV, where the Fjl^^j states 
are expected, our fits of the BaBar i?b-scan do not show a resonant structure due to the large decay widths of the 
states Y^hq]- resulting x^/d.o.f. = 88/67 with the 3 Breit-Wigners shown in Fig. [2] (right frame) is better than 
that of the BaBar fit. A Belle i?6-scan will greatly help to confirm or refute the existence of the state Yj;,^] visible 
in the analysis presented here. As the decays Y[f,(j] Bl Bg are not allowed, restricting the final states in to 
the Bq Bq (q — u,d), into which Y[;,g] decay, will reduce the background to the Yj^gj signal. It will be crucial to 
check that the characteristics of y[bq] (mass, full width and the electronic width) match those of the Fb, measured in 
the exclusive process e+e^ — > Yf, — > T(15, 25) tt+tt^. This may solve one of the outstanding mysteries in the T{nS) 
physics. 

The quantity TZee^Xb) in P3|) is given by the ratio of the two amplitudes Ay^^ ^ and Ayj^^^j, which also fixes the 
mixing angle 9. From our fit, we get 

7^ee(yb) = 1-07 ± 0.05, (38) 

yielding 

= -19±l° and AAf = 5.6 ± 2.8 MeV, (39) 

for the mixing angle and the mass difference between the eigenstates, respectively. 

The i?f,-analysis in the tetraquark picture can be used to determine k. The procedure how to do this requires 
some discussion, n can be determined from the theoretically estimated total decay widths of the Y[b^g] states and 
the corresponding result from the i?f,-fit. However, the estimated decay width of the Y[f,_g] is based on the input 
r[T(55)] = 110 ± 13 MeV from the PDG. The BaBar fit and ours, on the other hand, yield a lot smaller value for this 



14 



decay width (see, Table HX)) . To avoid the dependence on the absolute value of rtot[Tr(55')], it is safer to determine 
K from the ratios of the theoretical decay widths rtot(y[6,q])/rtot[Tr(5S')]t,hcory = (88 ± 16)k^/(110 ± 13), and the 
corresponding ratio of these widths obtained from the fit of the i?fc-data, Ttot{Y[b,q])/^tot[^{5S)]tit — (28 ±2)/ (46 ±8). 
This yields (adding the errors in quadrature): 



which is in the right ball park expected from the Lattice QCD estimates of the same [28|. For the mass eigenstates 
Y[b_i] and Y[h^h], the electronic widths ree{Y[t,j]) and ree{Y[b^h]) are given by Teeid) = 0.4 k'^Q{9)^ keV, as already 
stated. With the above determination of k and 9, we get 

Tee (Y[f,,;] ) = 0.09 ±0.03 keV and Tee (F[,,,^] )= 0.08 ± 0.03 keV. (41) 





11.1 11.2 
\fi [GeV] 



Y(10865) 



Yd 1020) 



V7 [OeV] 



FIG. 2: Measured Rb as a function of y/s with the result of the fit with 2 Breit-Wigners described in [ll (left frame). Reprinted 
from Fig. 1 of B. Aubert et al. [BaBar Collaboration], Phys. Rev. Lett. 102, 012001 (2009) [Copyright (2009) by the American 
Physical Society]. The result of the fit with 4 Breit-Wigners described in the text is shown in the right-hand frame, where we 
have indicated the location of the T(5S'), T(65') and the tetraquark state yjb,,] (labelled as F'^'). The location of the next 
higher J^'^ — 1 state Y^^]]] (labelled as Y'^') is also shown. The shaded bands around the mass of y(i' and y(2' refiect our 
theoretical uncertainty in the masses. 



TABLE IX: Fit values of the masses, decay widths (both in MeV) and the strong phases </> (in radians). 





M[MeV] 


T[MeV] 


ip [rad.] 


T(55') 


10864 ± 5 


46 ±8 


1.3 ±0.3 


T(65') 


11007 ±0.3 


40 ±2 


0.88 ± 0.06 




10900 - AM/2 ± 2 


28 ±2 


4.4 ±0.2 


Y\h,h] 


10900 4- AM/2 ± 2 


28 ±2 


1.9 ±0.2 



In conclusion, we have presented a case for the observation of a hidden bb tetraquark states in the BaBar Ri,- 
scan [l|. Our analysis is compatible with a J^^ — 1 state Yj^g] (10900) having a width of about 30 MeV. A scan 
of Rt in finer energy steps should be able to resolve the structure seen at this mass in terms of two mass eigenstates, 
split by about 6 MeV. The electronic widths are estimated to be between 50 and 120 electron volts. Other possible 
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371 and a dynamical model for the decays 



manifestations of tetraquarks have been discussed in the literature [36 
y6(10890) ^ T(15,25)7r+7r- is presented in |32|. 
Acknowledgments 

We thank the BaBar collaboration and the American Physical Society for their permission to show Fig. [2] (left 
frame) published in Helpful discussions with Riccardo Faccini and Antonello Polosa are gratefully acknowledged. 
We also thank Riccardo Faccini for providing us the fit program used in the BaBar analysis and we are grateful to 
Alexander Parkhomenko for reading the manuscript and pointing out several typos and notational inconsistencies in 
the first version of this paper. This work has been partially supported by funds provided by the ENSF, Trieste, Italy. 
Two of us (LA. and M.J. A.) would like to thank DESY for the hospitality during the summer 2009, where this work 
was done. C.H. wants to thank Benjamin Lutz for discussion on some technical aspects of the fits. 



[1] B. Aubert et al. [BaBar Collaboration], Phys. Rev. Lett. 102, 012001 (2009) [arXiv: 0809 .11201 [hep-ex]], and 

https:/ /oraweb. slac.stanford.edu/pls/slacquery/babar- documents. startup 
[2] K. F. Chen et al. [Belle Collaboration], Phys. Rev. Lett. 100, 112001 (2008) [arXiv:0710.2577] [hep-ex]]: L Adachi et al. 

[Belle Collaboration], arXiv:0808.2445 [hep-ex]. 
[3] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667, 1 (2008). 

[4] R. L. Jaffe, Phys. Rept. 409, 1 (2005) [Nucl. Phys. Proc. Suppl. 142, 343 (2005)] [arXiv:hep-ph/0409065| . 
[5] C. Quigg, Nucl. Phys. Proc. Suppl. 142, 87 (2005) [arXiv:hep-ph/0407124| . 
[6] A. Zupanc [for the Belle Collaboration]. larXiv:0910.3404l [hep-ex]. 

[7] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 74, 091103 (2006) [arXiv:he p-ex/0610018'. 
[8] M. Ablikim et al. [BBS Collaboration], Phys. Rev. Lett. 100, 102003 (2008) [arXiv:0712.1143 [hep-ex]]. 
[9] C. P. Shen et al. [Belle Collaboration], Phys. Rev. D 80, 031101 (2009) [arXl v:0808.0006 [hep-ex]]. 
[10] N. A. Tornqvist, Phys. Lett. B 590, 209 (2004) [arXiv:hep-ph/0402237| . E. Braaten and M. Kusunoki, Phys. Rev. 
D 69, 074005 (2004) [arXiv:hep-ph /031 1 147] . E. S. Swanson, Phys. Lett. B 588, 189 (2004) [arXiv7hep-ph/03ll229 



M. B. Voloshin, Phys. Lett. B 604, 69 (2004) [arXiv:hep-ph /040832l| . C. E. Thomas and F. E. Close, Phys. Rev. D 78, 

034007 (2008) ;arXiv:0805.3653 [hep-ph]]. 
[11] X. Liu, X. Q. Zeng and X. Q. Li, Phys. Rev. D 72, 054023 (2005) (arXiv:hep-ph/0507177 . X. Liu, Y. R. Liu, W. Z. Deng and 

S. L. Zhu, Phys. Rev. D 77, 034003 (2008) [arXiv:071 1 .04941 [hep-ph] ] . X. Liu, Z. G. Luo, Y. R. Liu and S. L. Zhu, Eur. Phys. 

J. C 61, 411 (2009) ;arXiv:0808.0073l [hep-ph]]. X. Liu and S. L. Zhu, Phys. Rev. D 80, 017502 (2009) [arXiv:0903.2529] . 
[12] J. L. Rosner, Phys. Rev. D 76, 114002 (2007) ,arXiv:0708.3496l [hep-ph]]. C. Meng and K. T. Chao, larXiv:0708.4222l 

[hep-ph]. S. H. Lee, A. Mihara, F. S. Navarra and M. Nielsen, Phys. Lett. B 661, 28 (2008) arXiv:0710.1029' [hep-ph]]. 

C. E. Thomas and F. E. Close, Phys. Rev. D 78, 034007 (2008) arXiv:0805. 3653 [hep-ph]]. N. Mahaian, arXiv: 0903.3107l 

[hep-ph]. T. Branz, T. Gutsche and V. E. Lyubovitskij, Phys. Rev. D 80, 054019 (2009) [arXiv:0903. 54241 [hep-ph]] . 
[13] E. Kou and O. Pene, Phys. Lett. B 631, 164 (2005) |arXiv:hep-ph/0507119i . F. E. Close and P. R. Page, Phys. Lett. B 



628, 215 (2005) [arXiv:hep-ph/0507199| . 



[14] L. Maiani, F. Piccinini, A. D. Polosa and V. Riquer, Phys. Rev. D 71, 014028 (2005) arXiv:hep-ph/0412098 . 

[15] L. Maiani, V. Riquer, F. Piccinini and A. D. Polosa, Phys. Rev. D 72, 031502 (2005) ,arXiv:hep-ph/ 0507062 . L. Maiani, 
F. Piccinini, A. D. Polosa and V. Riquer, AIP Conf. Proc. 814, 508 (2006) [arXiv:hep-ph/0512082| . L. Maiani, A. D. Polosa 
and V. Riquer, arXiv:0708.3997l [hep-ph]. L. Maiani, A. D. Polosa and V. Riquer, New J. Phys. 10, 073004 (2008). 

[16] N. V. Drenska, R. Faccini and A. D. Polosa, Phys. Lett. B 669, 160 (2008) |arXiv:0807."0593l [hep-ph]]. N. V. Drenska, 
R. Faccini and A. D. Polosa, Phys. Rev. D 79, 077502 (2009) arXiv:0902.2803J [hep-ph]]. 

[17] S. K. Choi et al. [Belle Collaboration], Phys. Rev. Lett. 91, 262001 (2003) arXiv:hep-ex/0309032l. 

[18] D. E. Acosta et al. [CDF II Collaboration], Phys. Rev. Lett. 93, 072001 (2004) arXiv:hep-ex/031202l] . 



16 



[19] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 93, 162002 (2004) [arXiv:hep-ex/0405004 



[20] B. Aubert et al. [BaBar Collaboration], Phys. Rev. D 71, 071103 (2005) arXiv:hep-ex/0406022 . 

[21] C. Bignamini, B. Grinstein, F. Piccinini, A. D. Polosa and C. Sabelli, Phys. Rev. Lett. 103, 162001 (2009) arXiv:0906.0882l 
[hep-ph]]. 

[22] P. Artoisenet and E. Braaten. [^Xiv:0911.20T6l [hep-ph]. 

[23] A. Abulencia et al. [CDF CoUaboration], Phys. Rev. Lett. 98, 132002 (2007) [arXiv:hep-ex/0612053l . 
[24] R. L. Jaffe, Phys. Rev. D 15, 281 (1977). R. L. Jaffe and F. E. Low, Phys. Rev. D 19, 2105 (1979). 
[25] R. L. Jaffe and F. Wilczek, Phys. Rev. Lett. 91, 232003 (2003) |arXiv:hep-ph/030734i; . 

[26] G. 't Hooft, G. Isidori, L. Maiani, A. D. Polosa and V. Riquer, Phys. Lett. B 662, 424 (2008) larXiv:0801 J288| [hep-ph]]. 

[27] M. G. Alford and R. L. Jaffe, Nucl. Phys. B 578, 367 (2000) [arXiv:hep-Iat/0001023) ; N. Mathur et al, Phys. Rev. D 
76, 114505 (2007) [arXiv:hep-ph/06071 lO] ; S. Prelovsek and D. Mohler, Phys. Rev. D 79, 014503 (2009) [arXiv:0810. 17591 
[hep-lat]]; H. Suganuma, K. Tsumura, N. Ishii and F. Okiharu, Prog. Theor. Phys. Suppl. 168, 168 (2007) [arXiv:0707.3309l 
[hep-lat]]. M. Loan, Z. H. Luo and Y. Y. Lam, Eur. Phys. J. C 57, 579 (2008) ;arXiv:0907.3609 [hep-lat]]; S. Prelovsek, 
T. Draper, C. B. Lang, M. Limmer, K. F. Liu, N. Mathur and D. Mohler, arXiv:0910.2749 [hep-lat]. 

[28] C. Alexandrou, Ph. de Forcrand and B. Lucini, Phys. Rev. Lett. 97, 222002 (2006) [arXiv:hep-lat/0609004| . 

[29] D. Ebert, R. N. Faustov and V. O. Calkin, Mod. Phys. Lett. A 24, 567 (2009) [arXiv:0812.3477l [hep-ph]] . 

[30] Z. G. Wang. larXiv:0908. 1266 [hep-ph]. 

[31] W. BuchmuUer and S. H. H. Tye, Phys. Rev. D 24, 132 (1981). 

[32] A. Ah, C. Hambrock and M. J. Aslam, DESY Report 09-222 [arXiv:0912.501 6' [hep-ph]]. 
[33] A. De Rujula, H. Georgi and S. L. Glashow, Phys. Rev. D 12, 147 (1975). 

[34] J. L. Domenech-Carret and M. A. Sanchis-Lozano, Comput. Phys. Commun. 180, 768 (2009) [arXiv:0805.2704l [hep-ph]]. 
[35] We thank Riccardo Faccini for providing us the of the BaBar fit shown in Fig. 1 of the BaBar paper ^ and the details 

of the modifications in replacing the flat nonresonant term by a threshold function at ^/s = 2ms . 
[36] W. S. Hou, Phys. Rev. D 74, 017504 (2006) ar Xiv:hep-ph/0606016| . 
[37] M. Karliner and H. J. Lipkin, arXiv:0802.0649, [hep-ph]. 



17 



